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Abstract — Space-time coding techniques have become 
common-place in wireless communication standards as they 
provide an effective way to mitigate the fading phenomena 
inherent in wireless channels. However, the use of Space-Time 
Block Codes (STBCs) increases significantly the optimal 
detection complexity at the receiver unless the low complexity 
decodability property is taken into consideration in the STBC 
design. In this letter we propose a new low-complexity decodable 
rate-1 full-diversity 4x4 STBC. We provide an analytical proof 
that the proposed code has the Non- Vanishing-Determinant 
(NVD) property, a property that can be exploited through the 
use of adaptive modulation which changes the transmission 
rate according to the wireless channel quality. We compare 
the proposed code to existing low-complexity decodable rate-1 
full-diversity 4x4 STBCs in terms of performance over 
quasi-static Rayleigh fading channels, detection complexity and 
Peak-to-Average Power Ratio (PAPR). Our code is found to 
provide the best performance and the smallest PAPR which is 
that of the used QAM constellation at the expense of a slight 
increase in detection complexity w.r.t. certain previous codes but 
this will only penalize the proposed code for high-order QAM 
constellations. 

Index Terms — STBCs, low-complexity decodable codes, condi- 
tional detection, nonvanishing determinants. 

I. Introduction 

SPACE-time coding techniques have become common- 
place in wireless communication standards (TJ, as they 
provide an effective way to mitigate the fading phenomena 
inherent in wireless channels. However, the use of Space- 
Time Block Codes (STBCs) increases significantly the optimal 
detection complexity at the receiver unless the low-complexity 
decodability property is taken into consideration in the STBC 
design. By complexity of detection, we mean the minimum 
number of times the receiver has to compute the Maximum 
Likelihood (ML) metric to estimate the transmitted symbols 
0. It is well known that Complex Orthogonal Design (COD) 
codes J4], are Single-Symbol Decodable (SSD) for general 
constellations. But if one considers the case of rectangular 
QAM constellations, the detection of each complex symbol 
reduces to separate detection of two real symbols which 
can be effectively performed via two threshold detectors (or 
equivalently PAM s) with a fixed complexity that does not 
depend on the constellation size. On the other hand, the 
rate of square COD codes decreases exponentially with the 
number of transmit antennas 0, which makes them more 
suitable for low-rate communications. In this letter we address 



the issue of increasing the rate of the COD at the expense 
of an increase in the detection complexity while preserving 
the coding gain and the Non- Vanishing-Determinant (NVD) 
property 0. This property can be exploited through the use 
of adaptive modulation which varies the transmission rate 
(through the choice of the modulation order) according to 
the wireless channel quality. The decoder exploits the STBC 
structure by estimating the set of the nx orthogonal symbols 
assuming the knowledge of the other set of n-i symbols, 
then performing the ML exhaustive search over symbols 
which reduces the complexity of detection from M^™ 1 " 1 "™ 2 ' 
to M" 2 where M is the rectangular QAM constellation size. 
Based on this conditional ML strategy 0, we construct a 
new low-complexity decodable rate-1 full-diversity 4x4 code 
that encloses the rate-3/4 COD code in 0. Furthermore, we 
provide an analytical proof (see Appendix) that the proposed 
code has the NVD property thus achieving the diversity- 
multiplexing trade-off 0. 

We compare the proposed code to existing low-complexity 
decodable rate-1 full-diversity STBCs in terms of performance 
over quasi-static Rayleigh fading channels, detection complex- 
ity and Peak-to-Average Power Ratio (PAPR). Our code is 
found to provide the best performance, the smallest PAPR, 
which is that of the used QAM constellation, at the expense of 
a slight increase in detection complexity w.r.t. certain previous 
codes only for large constellations. 

The letter is organized as follows: In Section II, preliminar- 
ies on STBCs are provided. In Section III we inttoduce the 
new code and derive its low-complexity decodability property. 
In Section IV, performance comparison by means of computer 
simulations of the proposed code and known similar codes in 
the literature is presented. We give our conclusions in section 
V, and the proof of the NVD property for the proposed code 
is provided in the Appendix. 

Notations: 

Hereafter, small letters , bold small letters and bold capital 
letters will designate scalars, vectors and matrices respectively. 
If A is a matrix, then A H , A T and det[A] denote the 
hermitian, the transpose and the determinant of A respectively. 
We define the vec(.) as the operator which, when applied to 
a m x n matrix, transforms it into a mn x 1 vector by simply 
concatenating vertically the columns of the corresponding 
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matrix. The eg) operator is the Rronecker product and 5kj 
is the kronecker delta. The sign(.) operator returns 1 if its 
scalar input is > and -1 otherwise. The round(.) operator 
rounds its argument to the nearest integer. The and 
operators denote the real and imaginary parts, respectively, of 
their argument. For any two integers a and b, a = b (mod n) 
means that a — b is a multiple of n. 

II. Preliminaries 
We define the MIMO channel input-output relation as: 

R = X H + W (1) 

TxN R TxN T N T xN R TxN R 

where T is the number of channel uses, Nr is the number 
of receive antennas, Nt is the number of transmit antennas, 
R is the received signal matrix, X is the code matrix, H is 
the channel matrix with entries hu ~ CAf(0, 1) and W is 
the noise matrix with entries w i3 ~ CJ\f(0,N o ). In the case 
of Linear Dispersion (LD) codes [9], an STBC that encodes 
2K real symbols is expressed as a linear combination of the 
transmitted symbols as: 

2K 

X = Y / PkX k (2) 

fc=l 

with Xk € R, K is the number of complex transmitted symbols 
and the /3/c, k = 1, 2K are T x Nt complex matrices called 
dispersion or weight matrices. For the LD codes, the MIMO 
channel model can be expressed in a useful manner as shown 
below. Expanding the code matrix X in the channel model as 
in (|2j, we obtain: 

2K 

R = ]T (fikH) Xk + W. (3) 

fc=i 

Applying the vec(.) operator to the R.H.S and the L.H.S we 
obtain: 

2K 

vec(R) = {In r ® fih) ™c (H) x h + vec(W). (4) 

k=l 

where Jjv b is the Nr x Nr identity matrix. 
If Yi, hi and w, designate the ith column of the received 
signal matrix R, the channel matrix H and the noise matrix 
W respectively, then equation Q can be written in matrix 
form as : 

r -i r ~i 
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Thus we have: 



r =Hs 



For COD codes, the weight matrices satisfy [5 |: 



(5) 
(6) 

(7) 



In this case, if is the i'th column of the equivalent channel 
matrix % in |6]), we may write: 

hf h, + hf h k = vec(H) H [(I Nr ® P k ) H (lN R ® Pi) 
+ (I Nr ® Pi) h {I Nr ® Pkj\ vec(H) 
= vec(H) H [I Nr ® MH) 
= 0, V jfc 7^ / 

where the last equality comes from Thus we have: 

3?{hfh/} = V k + I. (8) 

This relation can be used for the detection of orthogonal 
symbols as below. Expanding the MIMO channel model in 
(|6| we have: 



2 A' 



Us 



fc=i 



Multiplying both sides by we obtain: 



2K 



M\ 2 xi 



^hf h k x k + hf . 



(9) 



(10) 



fc=i 

k^l 



Taking the real parts of both sides, the second element of the 
R.H.S vanishes thanks to |8]l and we obtain: 

&{hfr} = ||vec(H)|| 2 x/ +5R{hfw} (11) 

■wi 

where wi is a zero-mean real Gaussian noise. If we assume that 
the complex symbols are drawn from a square constellation 
such as 4-/16-QAM constellations, which is usually the case 
in practical scenarios, the real symbols may be decoded 
independently through a PAM sheer as illustrated in (12) . It 
is worth noting that the PAM sheer equations require only 
a fixed number of simple arithmetic operations, which does 
not grow with the size of the square QAM constellation, and 
therefore they are considered of complexity 0(1). 

III. The Proposed Code 

From the previous section, a judicious structure for a new 
high-rate, low-complexity code would be to enclose a COD 
code into a new higher rate STBC. The resulting code will 
enjoy the low-complexity decodability through the use of 
conditional ML detection. 

The proposed code denoted X new takes the following form: 

X„ CH ,(s) = 0(an, ...,x 6 ) + (R 2 x 7 + R 3 x s ) R 4 (13) 
where the 0(xx, . . . , X2( a +i)) matrix denotes the square COD 
code for the case of 2 a transmit antennas and R^ is the 
matrix representations of the Clifford algebra generator 7^ 
with 1 < £ < 2a + 1, s = [si,S2], si = [x\,...,Xq\, 
s 2 = [^7)^8] an d 4> is chosen to maximize the coding 
gain. The code matrix can be explicitly written as in (T4| . 
It was found through exhaustive computer search that taking 
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Xi + JX 2 

-x 3 + jx 4 
-x 5 + jx 6 
-e^{-x 7 +jx 8 ) 



x 3 + jx4 x 5 + jx e -e^ixi + jxg) 

xi-jx 2 e^^xr+jxg) -x 5 -jx 6 

e^ixj + jx$) xi~jx 2 x 3 +jx 4 

X 5 - jx 6 -X 3 + jX4 xi + jx 2 



(14) 



4> = ^cos _1 (l/5) maximizes the coding gain and that it 
remains constant up to 64-QAM unnormalized constellations. 
An analytical proof that the code in ( fT4j > possesses the NVD 
property is provided in Appendix. The decoder exploits the 
structure of the proposed code by computing the ML estimates 
of the orthogonal symbols [x\, . . . , x§) assuming that a given 
value of {x*i : x%) is transmitted which we will denote s^ L |s2. 
To do so, the decoder scans all of the possible values of S2 
and assigns to each of these values a new signal z(s 2 ): 

8 

z(s 2 ) = r- ^2h k x k (15) 

k=7 

In the case of the square QAM constellations, the ML esti- 
mates of the orthogonal symbols (xx,...,Xg) assuming the 
knowledge of (x 7 , xg) can be obtained exactly as in ( p~2] > with 
the only difference that r will be replaced by z(s2). The 
decoder computes then the ML metric for the current value 
of (s^ L |§2,S2) according to: 



D(s 2 )= z(s 2 ) - ^ hi 



S2 



(16) 



The above metric calculation is repeated for every possible 
value of S2 which means that the decoder makes M metric 
computations. Finally, it decides in favor of the value of s 2 
with the minimum ML metric: 



_,ML 



arg mini) (§2) 



(17) 



where A is the set of the square QAM constellation points. 
The proposed decoder is illustrated in Fig. [T] 
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det and PAPR. The minimum determinant is defined as: 

Min det = min \det [(X(As))l I = y/8 (18) 

AseAC\{0} 

where As = s — s', AC is the vector space spanned by As and 
5 is the coding gain. For consistency, the minimun determinant 
is evaluated at fixed average transmitted power per channel use 
for all codes. The PAPR is defined as: 



PAPR„ 



max|X(i,n)| 2 
T-i£I£{|X(t,n)|2} 



(19) 



where t e {1,...T} and n e {1,...N T }. Due to the 
symmetry between transmit antennas, the subscript n will be 
omitted. It is worth noting that the complexity of detection 
of the QOD code in ifTUll can be reduced from 0{M 2 ) to 
O(M) and the complexity of detection of the Single Symbol 
Decodable (SSD) codes in lfTTI - lfT6l can be reduced from 
0(M) to 0(VM) as stated in flj) by the use of conditional 
detection and PAM sheers as for the proposed code. 

We can notice the following: The proposed code and the 
linear complexity decodable QOD code in ifTUl achieve the 
highest coding gain compared to existing low-complexity 
decodable rate-1 full-diversity 4x4 STBCs in the literature 
and thus we can expect that they provide better performance 
in the high SNR region which can be verified from Figs. [2] [3] 
The proposed code along with the codes of [10], ifTTI have the 
lowest PAPR which is that of the used constellation. On the 
other hand, our code suffers from a slight increase in decoding 
complexity w.r.t. the codes of ifTTl - lfTBI . but this will only 
penalize the proposed code for high-order QAM constellations 
(in fact our code will be more complex to decode only for 
M > 64). The performance of the proposed code in terms of 
Codeword Error Rate (CER) versus SNR per receive antenna 
is compared with other low-complexity decodable rate-1 full- 
diversity STBCs in the literature. Simulations are carried out 
in a quasi-static Rayleigh fading channel in the presence of 
AWGN for 1 and 2 receive antennas and 2 bits per channel use 
(bpcu) (see Figs. 2|3 1. One can easily verify that the proposed 
code and the Quasi-Orthogonal Design (QOD) code in IfTUl 
provide the best performance. The rate-1 4x4 SAST, EOS 
and EAST codes in lfT3l . IT31 . |[T6l respectively, provide the 
same performance as the STBCs in [12|, [14| and therfore 
have been omitted to avoid congestion in Figs. [2] [3] 



Fig. 1 . The low-complexity decoder of the proposed code 

IV. Numerical and Simulations Results 

Table [I] summarizes the comparison of the proposed code 
with existing similar low-complexity decodable codes in terms 
of detection complexity for square QAM constellations, Min 



V. Conclusion 

In the present letter, we have introduced a new low- 
complexity decodable rate-1 full -diversity 4x4 STBC that 
encloses the rate-3/4 COD and retains its coding gain. The 
coding gain is proved analytically to be constant over QAM 
constellations. Simulation results have shown that the perfor- 
mance of the proposed code is identical to the performance of 
the code in iflOll while being decodable at a lower complexity 
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TABLE I 

SUMMARY OF COMPARISON IN TERMS OF COMPLEXITY, MlN det AND PAPR 



Code 


Detection 
Complexity 


Min det (= \/S) 
for QAM constellations 


PAPR (dB) 


QPSK 


16QAM 


64QAM 


The proposed code 


M 


16 





2.55 


3.68 


The QOD code in [10] 


2M 


16 





2.55 


3.68 


The SSD code in HU 




7.11 





2.55 


3.68 


The codes in [12], [13] 


aVm 


12.8 


1.61 


4.16 


5.29 


The codes in [14|, [15|, [16] 


4VM 


12.8 


5.79 


8.34 


9.47 
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— * — Proposed code 
-V — The code in [10] 
-6 — The code in [1 1 
-A — The code in [12] 
— b — The code in [1 4] 
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Fig. 2. CER performance for 4x1 configuration and 2 bpcu 
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Fig. 3. CER performance for 4x2 configuration and 2 bpcu 

and it outperforms the codes in lfTTI - lfT6l at the expense of 
a slight increase in decoding complexity for the 64-QAM 
constellation and above. The proposed code has the same 
PAPR as the QOD code in (I0l and the code in ifTTl . 

Appendix 

In the following we will prove that choosing = 
| cos -1 (1/5) indeed maximizes the coding gain 6 (which 
is equal to 16 2 ) for unnormalized QAM constellations and 



guarantees the NVD property for the proposed code. The 
coding gain 6 is equal to the minimum Coding Gain Distance 
(CGD) HI, or mathematically: 



S = min det 

s.s'eC 



(X(s)-X(s')) ff (X(s)-X(s')) 



CGD(X(s),X(s')) 



min \det[(X( As))} \ 2 

AseAC\{0} 



(20) 



where As = s — s', AC is the vector space spanned by As. 
The code structure in ( fT4| > imposes: 

min \det [(X(As))l I < \det [0(2, 0, 0, 0, 0, 0] I = 16 

AseAC\{0} 

As a result, the angle <f> that maximizes the coding gain has 
to satisfy: 

\det [(X(As))] | > 16, V As ^ (21) 
For the proposed code we have: 

/ 6 \ 2 / 8 \ 2 

|<fcr[(X(As))]| 



^Ax 2 +e^b+e^ [^Aa 



(22) 



where: 



+ (Ax 7 Ax 6 + Ax 8 Ax 2 ) 2 + (Ax 7 Ax 4 - Ax 8 Ax 3 ) 2 . 

A simplification of the expression of b is possible by noting 
that: 

(E Aa; *) (E Ax * ) - ( Ax 2 + Ax i) (A^ 2 7 + Ax 2 ) 

+ (Ax 2 + Ax 2 ) (Ax 2 + Ax|) + (Ax 2 + Ax 2 ) (Ax 2 + Ax|) 
Applying the Fibonacci's two square identities we obtain: 

(Y1 Ax ^J i^2 Ax2 }j = (Ax 7 Ax 2 - Ax 8 Ax 6 ) 2 

+(Ax 7 Ax 6 + Ax 8 Ax 2 ) 2 + (Ax 7 Ax 4 - Ax 8 Ax 3 ) 2 
+(Ax 7 Ax 3 + Ax 8 Ax 4 ) 2 + (Ax 8 Axi - Ax 7 Ax 5 ) 2 

6 

+(Ax 7 Axi + Ax 8 Ax 5 ) 2 = J^a 2 , a, = 45i, a,* G Z. 

»=i 

Therefore, one may write b in more compact form a: 



(23) 



5 



Setting x = e) 2 ^, the discriminant A of the second degree 
equation |22) is expressed as: 




) 4 



A = 4 



Consequently, the roots of equation ( |22| ) are: 




Al,: 



b 2 



Ax: 



(25) 



where A 2 = A*, |Ai| = |A 2 



>:; - v 

For the sake of simplicity, we will denote hereafter 
Y^i=\ Ax 2 = a i and 52j- =7 Ax 2 = cr 2 . In the case of <r 1 ^ er 2 , 
equation. [22] can be lower bounded as below: 



\clet [(X(As))] 



As^O 



CI 



(x - Ai) (a; - A 2 ) 



><t| 



(M 



*2 



cr 2 



N)(W 

2 



|A 2 |) 



(cr 2 - (Ti) 



> 16 
(26) 

where the latter inequality follows by substituting Axi — 
2n i: rii E Z as we are dealing with unnormalized QAM 
constellations. If o\ = cj 2 = a, equation ( |22] > can be written 

as 



|^[(X(As))]| 

where b = 2 (cr 2 - 2(a? 
we have to prove that: 

2a 2 



As^O 



= |2ct 4 cos(20) + b\ (27) 
a|)). Taking cos(20) = 1/5, 
(28) 



> 16 V ai = <t 2 = a 7^ 0. 

Multiplying both sides by 5 and using ((23), the above inequal 
ity becomes: 



Il2cr 2 - 20(a 2 



*2Tu. 3 )|( ai ,..., Q6)7 io 

Denoting a* = ^ and cr = ? we may have: 

|i2a 2 ~20(a 2 + a 2 + ^)| ((ii _ 

The above inequality is satisfied iff: 

|3a 2 -5(a 2 + a 2 + a 2 ) l(5i rae) ^ Q 

However, the L.H.S of the above inequality can be considered 
as a special case of: 

3X 2 - 5 (X 2 + X 2 +X 2 ) (32) 

where X, t E Z. This type of equations have been extensively 
studied in the mathematical literature dealing with the solv- 
ability of quadratic Diophantine equations (see |[T9l ). Applying 
theorem 6 in lfl9ll we have: 

3X 2 - 5{Xl+X 2 +Xl) + V (X 1 ,X 2 , X 3 ,X 4 ) ? (33) 



4) 



> 5 x 16 



(29) 



> 5. 



> 2. 



(30) 
(31) 



as —3 x 5 x 5 x 5 
4 (mod 8). 
Moreover, 



1 (mod 8) with 3-5-5-5 



3X 2 - 5pf 2 + X 2 + Xf) ± ±1. 



-12 



(34) 



Otherwise, we must have 

-2 



3X'f = ±1 (mod 5) (35) 

which cannot be true, since the quadratic residues modulo 
5 are 0,1 and 4 ED, thus 3X 2 = 0, ±3 or ± 2 (mod 5). 
Therefore, we can write: 



3X 2 - 5{X 2 + X 2 + Xj) 
which in turns implies: 

|<fct[(X(As))] 



(X 1 ,X 2 ,X 3 ,X i ) 9 iO 



> 2 



Ast^O 



> 16 



(36) 



(37) 



thus ending the proof. 
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